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43 ■ We investigate the phase diagram of the so-called Polyakov-Nambu- 

Jona-Lasinio (PNJL) model at finite temperature and nonzero chemical 
potential. The calculations are performed in the light and strange quark 
sectors (u, d, s), which includes the 't Hooft instanton induced interaction 
term that breaks the axial symmetry, and the quarks are coupled to the 
(spatially constant) temporal background gauge field. On one hand, a spe- 
cial attention is payed to the critical end point (CEP). The strength of the 
flavor-mixing interaction alters the CEP location, since when it becomes 
weaker the CEP moves to low temperatures and can even disappear. On 
the other hand, we also explore the connection between QCD, a nonlocal 
Nambu-Jona-Lasinio type model and the Landau gauge gluon propaga- 
tor. Possible links between the quenched gluon propagator and low energy 
hadronic phenomenology are investigated. 

PACS numbers: ll.10.Wx, 11.30.Rd, 12.40.-y 

1. The phase diagram in the context of the PNJL model 

Chiral symmetry breaking and confinement are two of the most impor- 
tant features of quantum chromodynamics (QCD). Chiral models like the 
Polyakov-Nambu-Jona-Lasinio (PNJL) model have been successful in ex- 
plaining the dynamics of spontaneous breaking of chiral symmetry and its 
restoration at high temperatures and densities/chemical potentials. The 
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PNJL model also plays an interesting role in the investigation of the QCD 
phase structure. Understanding the properties of matter at finite tempera- 
tures and densities is one of the most important goals from both the theo- 
retical and experimental point of view. For example, the critical end point 
of QCD, proposed at the end of the eighties, is still a very important subject 
of discussion nowadays: indeed its existence and location is one of the main 
goals in SPS at CERN and in RHIC at BNL [I]. 

The NJL model describes interactions between constituent quarks, giv- 
ing the correct chiral properties; static gluonic degrees of freedom are then 
introduced in the NJL lagrangian, through an effective gluon potential in 
terms of Polyakov loops, with the aim of taking into account features of both 
chiral symmetry breaking and deconfinement. The coupling of the quarks 
to the Polyakov loop leads to the reduction of the weight of quark degrees 
of freedom as the critical temperature is approached from above, which is 
interpreted as a manifestation of confinement and is essential to reproduce 
lattice results. 

Our calculations are performed in the framework of an extended SU(3) f 
PNJL Lagrangian, which includes the 't Hooft instanton induced interaction 
term that breaks the U^(l) symmetry, and the quarks are coupled to the 
(spatially constant) temporal background gauge field $ [21 [3]: 

1 8 

C = q(i^D tM -m)q+-g S £ [ ( q X a q f + ( qi 75 X a q f ] 

1 a=0 

+ g D {det[q(l + 75 )q] + det[q(l- l5 )q}}-U($[A},$[A}-,T). (1) 

The covariant derivative is defined as = — iA^, with A^ = 5qAq 
(Polyakov gauge); in Euclidean notation Aq = —1A4. The strong coupling 
constant g is absorbed in the definition of A fJ, (x) = gA^(x)^ , where is 
the (SU(3) C ) gauge field and A a are the (color) Gell-Mann matrices. 

The effective potential for the (complex) field adopted in our parame- 
trization of the PNJL model reads: 
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The parameters of the effective potential U are given by ao = 3.51, 
ai = —2.47, ci2 = 15.2 and ^3 = —1.75. When quarks are added, the pa- 
rameter Tq, the critical temperature for the deconfinement phase transition 
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(that manifests itself as a breaking of the center symmetry) within a pure 
gauge approach, was fixed to 270 MeV, according to lattice findings. This 
choice ensures an almost exact coincidence between chiral crossover and de- 
confinement at zero chemical potential, as observed in lattice calculations. 

The parameters of the NJL sector are: m u = = 5.5 MeV, m s = 140.7 
MeV, g s A 2 = 3.67, g D A 5 = -12.36 and A = 602.3 MeV, which are fixed to 
reproduce the values of the coupling constant of the pion, f n = 92.4 MeV, 
and the masses of the pion, the kaon, the rj and rf, respectively, M n = 135 
MeV, M K = 497.7 MeV, M v = 514.8 MeV and = 960.8 MeV @|. 

The inclusion of the Polyakov loop effective potential Z//(<5, <J>; T), that 
can be seen as an effective pressure term mimicking the gluonic degrees 
of freedom of QCD, is required to get the correct Stefan-Boltzmann limit. 
Indeed in the NJL model the ideal gas limit is far to be reached due to the 
lack of gluonic degrees of freedom. 

In Fig. [1] (left panel), we present the phase diagram of the PNJL model. 
As the temperature increases the chiral transition is first order and persists 
up to the CEP. At the CEP the chiral transition becomes a second order 
one. The location of the CEP is found at T CEP = 155.80 MeV and fi CEP = 
290.67 MeV {p% EP = 1.87p )- For temperatures above the CEP there is a 
crossover whose location is calculated making use of d 2 (qq) /dT 2 = 0, i.e. 
the inflection point of the quark condensate (qq) . 

The transition to the deconfinement is given by d 2 Q/dT 2 = 0, and is 
represented by the magenta line. The surrounding shaded area that limits 
the region where the crossover takes place is determined as the inflection 
point of the susceptibility dQ/dT. 

Due to the importance of the location of the CEP from the experimental 
point of view, let us investigate the influence of other parameters which can 
lead to a significant change in the CEP's localization. 

It is well known that the Ua(1) anomaly has big influence on the behav- 
ior of several observables, so it is demanding to investigate possible changes 
in the location of the CEP in the (T, /i) plane when the anomaly strength 
is modified. The axial Ua(1) symmetry is broken explicitly by instantons, 
leaving a SU(N/)(g> SU(Nj) symmetry which determines the chiral dynam- 
ics. Since instantons are screened in a hot or dense environment, the Ua(1) 
symmetry may be effectively restored in matter. So, the change of the Ua(1) 
anomaly strength has a strong influence on the localization of the CEP in 
the (T, fx) plane. 

In Fig. [1] (right panel), we show the location of the CEP for several 
values of go compared to the results for gjj , the value used for the vacuum. 
As already pointed out by K. Fukushima in j2j, the location of the CEP 
depends on the value of go- The results show that, in the framework of this 
model, the existence or not of the CEP is determined by the strength of 
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Fig. 1. Phase diagram in the SU(3) PNJL model: the location of the CEP is found 
at T CEP = 155.80 MeV and fi CEP = 290.67 MeV (see details in text). Dependence 
of the location of the CEP on the strength of the 't Hooft coupling constant go- 



the anomaly coupling, the CEP getting closer to the \x axis as go decreases. 
As the strength of the flavor-mixing interaction becomes weaker, the CEP 
moves to low temperatures and can even disappear. 



2. Low energy physics and the gluon propagator 

In this section we explore the connection between QCD, a nonlocal 
Nambu-Jona-Lasinio type model and the Landau gauge gluon propagator 

m. 

The interaction between quarks and gluons in QCD reads: 

C&A = rf-fA%^* (4) 

Expanding the term containing C^ A up to g 2 and integrating the gluon 
fields (see [5] for details), the theory becomes an effective nonlocal fermionic 
theory 

S[tp, ip] = J d 4 xd 4 y{ip(y) 5{y - x) (ry^ - m) ip(x) 

+ -^J{x,y)D{x - y)J{y,x) - — J 5 (x,y)D(x - y)Js(y,x)j (5) 

with J(x, y) = tp(x)ijj(y) and J$(x, y) = tp(x)j^ip(y) and D(x, y) is the gluon 
propagator form factor. 
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First principles calculations of the gluon propagator have been per- 
formed using lattice QCD and DSE (see for example [6] and references 
therein). The momentum space propagator 

D(p 2 ) = Z^n (6 ) 

Was 



\r + aqcd 

is able to describe both the scaling (re > 0.5) and decoupling (k = 0.5) 
infrared DSE solutions and the lattice data up to p ~ 800 MeV; Aqcd 
stands for an infrared mass scale. 

Let us define the dimensionless form factor in momentum space as 




iV) = A 2 D(p 2 ) = - — 9(A-p). (7) 



The constant Z in Eq. ([6]) will be included in the definition of the coupling 
constant G, which multiplies the quark currents of the nonlocal theory. G 
carries the dimension of a length squared. In f(p 2 ), A is the cut-off. In a 
first step we assume Aqcd = A. The form factor f(p 2 ) is shown in Fig. [2] 
(left panel), together with typical form factors considered in the literature. 

Demanding the nonlocal model to reproduce the experimental values 
for M n , and IV-vy^, with a cut-off A = Aqcd = 800 MeV, we obtain 

m q = 4.205 MeV, - (qq) 1/S = 271.1 MeV, GA 2 = 7.491 and k = 0.529. The 
presented results favor k > 0.5. Now we investigate the decoupling type of 
propagator 

Z 



Luon 



where M g i uon takes the role of Aqcd and can be interpreted as an effective 
gluon mass. Requiring the model to reproduce the same experimental quan- 
tities as before, we find M g i uon within typical values found in the literature, 
but with a strong dependence with the cut-off. In fact, M g i uan is a linear 
function of A - see Fig. [2]- right panel. In what concerns the quark conden- 
sate, the model shows that (qq) increases with M g i uon , i.e with the cut-off A. 
In order to reproduce the experimental value of the condensate, i.e. to have 
(- (<?<?)) 1/3 = 270 MeV, it turns out that the gluon mass is M g i uon — 878 
MeV for A = 800 MeV. We therefore conclude that low energy physics does 
not distinguish between the so-called decoupling and scaling solutions of 
the Dyson- Schwinger equations. This result means that, provided that the 
model parameters are chosen appropriately, one is free to choose any of the 
above scenarios. 

Finally, it is interesting to refer that the model considered here is chiral 
invariant and satisfies the GMOR relation at the 1% level of precision. 
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Fig. 2. Form factors as a function of p. The figure includes typical form factors used 
in previous studies. M g i uon (A) required to reproduce the experimental rV-^p,. 



As future work we will generalize the results to non-zero temperature 
(this requires modeling the gluon propagator at finite temperature by a 
functional form compatible with both Dyson-Schwinger and lattice QCD 
results) which will allow us to investigate the meson properties at finite 
temperature as probes for the chiral symmetry restoration. 
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